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ABSTRACT 

A new method for detecting negative cycles in a graph is 
proposed. This method is based upon the primal - dual 
relationships of a linear program formulated from an assign- 
ment problem type network. A computer program is developed 
for this new method to include the complete solution of the 
assignment problem. Results are given on program efficiency. 
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I. OBJECTIVE 



A new method has been proposed by H. Greenberg for the 
detection of negative cycles in an assignment type problem 
The objective of this thesis was to develop a computer pro 
gram for this new method and to measure the efficiency 
thereof relative to execution time. 
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II. DESCRIPTION OF THE METHOD 



A. INTRODUCTION 

Suppose there exists a network with vertices V = {l,...,n}, 
directed edges E = {(i,j)e V X V}, a non-negative integral 
valued function k giving the maximum allowable flow k— over 
every edge, and a non-negative cost function a giving the cost 
a^j associated with a unit of flow over any edge. A flow X 
of value v is an integral valued function on E satisfying: 



( 1 ) 0 < x. - < k. . , (i, j)e E 

K J _ !-)_ 13 ’ v * JJ 



(x. . - x..) = v, i = 1 , 

= 0 , i = 2 , n- 1 , 

= -v, i = n. 

Vertices 1 and n are, respectively, called the flow source 
and sink. 

The minimal cost flow problem is to find, among all flows 
X of value v, one which minimizes: 

(3) Q(X) = E x. • a. . . 

(i,j)e E J 3 

The assignment problem is to fill n jobs by as many men 
at least total cost where a^j represents the cost of using 
man i (i = 1, . . . , n) in job 3 (j = n + 1, ..., 2n) . Because 
of this bipartite form the problem can be thought of as a 
special minimal cost flow problem. 

Now suppose a flow has been found satisfying (1) and (2) , 
and for the assignment problem this corresponds to finding a 
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Rule" given by Dantzig in [1] (p. 307). We then construct an 
associated network G(X), given by Klein [2], which has the 
same vertices as the original network and directed edges as 
follows : 



These revised edge costs are simply those associated with 
increasing or cancelling the flow by one unit on these edges: 
the revised capacities indicate the extent to which this can 
be done. 

Now, use a result proved in Busacker and Saaty ([3], pp. 
256-257) . 

Theorem . X is a minimal cost flow if and only i:f there 
is no directed cycle C in G(X) such that the sum of the costs 
around C's edges are negative. (A directed cycle is a 
sequence of distinct directed edges of the form 
{(io’ii) ( i i» i 2 ) ' ** Ci p ,i q ) (i q >i 0 )} involving distinct 
vertices . ) 



(4) E(X) : (i , j) if x^ < k^ and x^ = 0 
(j ,i) if >0, 



with revised capacities: 





. . = x. . , 
il il 




and with revised edge costs: 



(6) a.' : a.". = a. . , if x. . < k. . and x.. = 0 

il il il il li 
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An immediate consequence of this theorem is that a test 
of the optimality of the flow X is at hand if G(X) can be 
checked for the existence of a negative cost directed cycle. 
Further, if such a cycle is found, an improved flow is ob- 
tained by sending a positive unit flow around this cycle. 
Such a flow alteration obviously leads to a lower total cost 
and also leaves the flow value v unchanged. 



B. THE NEW METHOD 

Several methods exist for locating negative cost directed 
cycles (see Klein [2]). The new method adds to these. The 
idea behind this method is to formulate a linear program 
from the network G(X) using the flow balance equations as 
constraints. A mathematical statement of this set of 
equations for the assignment problem is: 



(7) minimize: Z 



£ x. . 
(i,j)eE 1J 




where E(X) is deter- 
mined by the assignment, 



subject to: 



E x. . - E *ij = 0 , i=l, • • • ,n, 
j > ( i > j ) e E 3 j,(j,i)eE 

Flow J 

Balance ) E x. . - E x ij = 0 , j =n+l , • • • , 2n , 

l i > C j > i) £ E 13 i,(i,j)eE 



Requires 

positive 

flow 



n 2n 

E E X- • 1 , 

i=l j=n+l 1J 



Obviously, a feasible solution exists to this program. 
But, more importantly, if a negative circuit exists then the 
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solution will be unbounded since the sum of the a around 
the circuit will be negative and the flow around this circuit 
will be increased with the solution approaching minus infinity. 

The consequence of this is that we have at hand another 
method for detecting negative circuits by capitalizing on the 
theorem (see Hadley [4]) that an infeasible dual implies an 
unbounded primal when the primal has at least one feasible 
solution. Thus the dual of the above program (7) is solved 
through Phase I, only, of the Simplex Method to test for 
feasibility. If feasible, then the assignment is optimal. 

If the dual is infeasible, then the assignment is not optimal 
and a change of assignment is made. This process is repeated 
until a feasible dual and, hence, an optimal assignment is 
obtained . 

The formulation of the primal program is now illustrated 
with a two by two assignment problem example. The cost matrix 
for this example is (8) : 

( 8 ) 

3 4 

1 
2 

where i = 1,2 denote the men and j =3,4 denote the jobs, 
hereafter referred to as assignee and assigned nodes, 
respectively. First, an initial assignment must be obtained, 
say by the "Least-Cost Rule". Suppose this assignment is 



a 13 


a 14 


a 23 


a 24 
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(j-i), 3-1, 4-2. Then the graph G(X) representing this 
problem can be constructed (9) . 

(9) 



" a 13 X 13 




The flow balance equations are then (10) : 

( 10 ) 



Node 


1: 


X 14 


X 13 > 


Node 


2: 


X 23 


X 24 ’ 


Node 


3: 


X 13 


X 23’ 


Node 


4: 


X 24 


= x 14 * 



From which is written the linear program (11) : 

(11) minimize: Z = ‘ a 1 3 x 13 + a i 4 x i 4 + a 23 x 23 " a 24 X 24’ 

subject to: 



X 13 


+ X 14 








o. 








X 23 


X 24 


= 


0 


> Flow 


X 13 




~ X 23 




= 


0 , 


j Balance 




X 14 




+ X 24 


zz 


°> j 


1 


X 13 


+ X 14 


+X 23 


+ X 24 


> 


1. 1 


l Positive 




x. . 
il 


> 0 . 






J 


Flow 
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The dual to this program (11) is then solved through Phase I 
to test for optimality of the assignment. 
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III. PROGRAMMING THE ALGORITHM 



A. 



THE STANDARD FORMAT 

The dual program to (11) is (12) : 

(12) maximize v^, 



subject to: 



-v. 



v. 



-v. 



+v. 



■ V- 



- a 






+v y 



< 

< 

< 



13’ 
l 14 ’ 

1 23 ’ 

1 24 ’ 



which, after conversion to the standard form for solution 
is (13) : 

(13) maximize , 



subject to: 



v. 



V, 



-v. 



V 2 " V 3 



= a 



-v, 



13’ 

14 



' s. 



V, 



- V, 



3 23 

2 v 4 ” s 4 = a 24 ’ 

with slack variables s^ through s^. Examining the array of 
coefficients, a standard format for any size problem can be 
detected. This format can be illustrated with a three by 
three problem annotated with general dimensions (14) , where 
n denotes the number of assignee or assigned nodes (in this 

2 

case n = 3), v^ (i-l , • • • , 2n) the dual variables, Sj(j=l, , **,n ) 
the slack variables, -s^ the artificial variables added to 
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( 14 ) 







-6rn 


<S-n 










ii >> 






Basic 

Variables 


V 1 V 2 V 3 


V 4 V 5 V 6 


S 1 S 2 S 3 S 4 S 5 S 6 S 7 S 8 S 9 


Constants 


A 


S 1 


1 


-1 


1 


a 14 


n 


" S 2 


1 


-1 


-1 


a 15 


I 


S 3 


1 


-1 


1 


a 16 




~ S 4 


1 


-1 


-1 


a 24 


2 1 
n n 


S 5 


1 


-1 


1 


a 25 




s 6 


1 


-1 


1 


a 26 


t 


S 7 


1 


-1 


1 


a 34 


n 

i 


S 8 


1 


-1 


1 


a 35 


y v 


" S 9 


1 


-1 


-1 


a 36 




-w 


-1-1-1 


111 


010 1 00001 


- (a 15 +a 24 +a 36 ) 



obtain a basis, and W the infeasibility form. Examining the 
format in detail, it is noted that the first n columns derive 
from the flow balance constraints on the assignee nodes 
(numbered from 1 to n) ; the second n columns derive from the 
constraints on the assigned nodes (numbered from n+1 to 2n) ; 
and the n slack columns reflect the assignment. The slack 
columns reflect the assignment in that the rows in which the 
-l’s appear indicate the assignment that has been made. If 
the rows in each group of n rows are numbered identically from 
n+1 to 2n then the row number in which the -1 appears in the 
first group of n rows is the number of the node assigned to 
the first assignee node. Similarly, the remainder of the 
assignment can be determined. In this example (14) , 5 is 
assigned to 1, 4 is assigned to 2, and 6 is assigned to 3. 



12 



Of course, only n of these -l's may appear. Moreover, the row 
number of a -1 may not be duplicated between groups since this 
would indicate a duplicate assignment. 

It is seen then that the only format change between dif- 
ferent problems is in the slack columns reflecting the 
assignment. This feature was exploited in programming the 

simplex iterations by reducing the computer storage require- 

2 

ment for this matrix to a single vector of dimensions n by 
1. This vector records the assignment and its contents are 
scanned when format information is needed. It should also be 
noted that with this format the matrix of coefficients -is 
unimodular. Hence, for every iteration the coefficients 
will be +_ 1 or 0 . 

B. ITERATIVE PROCEDURE USED 

The revised simplex method was used to perform the Phase 
I iterations because of the computational advantages. Exam- 
ining the standard format (14) , we see that at least n 
iterations are required because of the -l's in the first n 
rows of the infeasibility form equation. To reduce compu- 
tational time the revised simplex method was modified to carry 
out these n iterations. The simplex multipliers were not used 
to find the candidate for the new basic variable. Instead, the 
first n variables were introduced into the basis sequentially. 
Subsequently, the standard revised simplex method was used. 

At this point it would be useful to illustrate the contents 
of the matrix used to perform the iterations. Again, this will 
be illustrated with the previous three by three problem, anno- 
tated with general dimensions (15) . This matrix is denoted B 
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( 15 ) 






n 



A 



n 



S l’ S 2 S 3" S 4 S 5 S 6 S 7 S 8’ S 9 



Constants 



14 



15 



16 



24 



25 



26 



34 



35 

L 36 



ACOL 

0 

1 

0 

1 

0 

0 

0 

0 

1 



VCOL 

7 

-8 

9 

-10 

11 

12 

13 

14 
-15 



0 - 10-1000 0 -1 



" ( a 15 +a 24 +a 34^ 



in the computer program. In addition to containing the in- 
verse of the basis and the simplex multipliers, where - s ^ 
denotes the artificial variable added for the corresponding 
negative slack variable, the matrix contains an assignment 
column (ACOL) and a variable identification column (VCOL). 

The assignment column, with binary entries of 0 or 1 only, 
where l’s record the assignment (in this case the initial 
assignment is 1-5, 2-4, 3-6), is used as a memory device in 
performing the iterations and in printing the solution. The 
variable column records the variables in the basis and is up- 
dated after each iteration. The minus indicates an artificial 
variable. If a problem is infeasible, this column is scanned 
to determine which artificial variables remain in the basis. 
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This information is necessary in changing the assignment, 
which shall be described later. 

While the revised simplex method requires less computer 
storage than the standard simplex method, the storage require- 
ments for this method are, nonetheless, relatively large. For 
instance a 15 by 15 problem requires approximately 210,000 
bytes of storage. (228 X 227 = 51,756 storage locations for 
the matrix B + 675 storage locations for a computing register 
and input/output requirements = 52,431 total X 4 bytes per 
location for single precision = 209,724 bytes.) But, other 
methods for detecting negative circuits require as much 
storage . 

C. INITIAL ASSIGNMENT 

Up to this point an initial assignment has been assumed, 
but the method for finding that assignment has not been 
described. The method selected was the "Least-Cost Rule" 
with steps as follows : 

1) The cost matrix is scanned to determine the lowest 
cost, i.e. , find Min a. . . 

(i,j) 13 

2) The column j and the row i in which this cost 
appears constitutes the first assignment. That column and 
row is then ineligible for further assignment. 

3) Steps 1) and 2) are repeated for costs remaining 
in eligible rows and columns until all assignments have been 
made . 
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D. CHANGE OF ASSIGNMENT RULE 



Summarizing the procedure described thus far - First, an 
initial assignment is determined. Next the dual problem to 
(7) is solved through Phase I of the simplex method. If the 
problem is feasible, then the initial assignment is optimal 
and the solution is printed. If the problem is not feasible 
then the initial assignment is not optimal and a change in 
assignment is required. 

The key to devising an assignment change rule is contained 

in the artificial vectors remaining in the basis because they 

2 

reflect the infeasibility. In every column of the last n 
columns in which an entry other than zero appears in the 
artificial vector, the sign is changed of the + 1 in that 
column of the basic format (14). Thus, the artificial vectors 
remaining in the basis are identified and constructed from the 
inverse and the procedure just described is performed for each 
of these vectors sequentially. The efficacy of this rule is 
evaluated in the following section. It should be noted that 
since this procedure does not change the coefficients of any 
variable in the basis, the inverse for the last tableau remains 
valid. Therefore, it is not necessary to cycle back to the 
initial tableau. Rather, after the assignment has been recorded 
in the memory device (ACOL) , the simplex iterations can continue 
from the last tableau. 

A good illustration of application of the assignment rule 
is problem 5 (see Appendix B) , a five by five problem. After 
nine iterations the program was found to be infeasible and two 
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artificial vectors remained in the basis (-s^g and -s^) . 

The vectors shown on the following page trace the change in 
assignment from initial contents of ACOL to final contents 
of ACOL. For ease in recognizing the assignments, the vectors 
have been partitioned by n (n = 5) . Recall that the entries 
of ACOL are binary (0,1) and that an entry other than 0 in 
the artificial vector "flip-flops" the corresponding ACOL 
entry . 
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IV. RESULTS 



Twenty problems, ranging in size from four by four to 
15 by 15, were tried. The entires comprising the cost matrix 
were the integers one through nine, selected randomly (see 
Appendix B) . The method yielded an optimal solution in each 
case. The execution times ranged from 0.7 seconds for a 
four by four to 56 seconds for a 15 by 15 (see Appendix A) . 

In three instances the initial solution was optimal and, 
hence, an assignment change was not required. But, in the 
remaining 17 problems an assignment change was required and 
in every case the change rule produced the optimal assignment 
on the first attempt. With respect to the number of iter- 
ations required to detect a negative circuit the range was 
from 7 to 33. And, in all cases the number of iterations 
required was less than the number equations involved. 
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APPENDIX A -TIMING RESULTS 
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*Initial problem was feasible; i.e., initial assignment was optimal 



APPENDIX B 



PROBLEMS AND SOLUTIONS 



1. Problem 1 

Initial Assignment Optimal Assignment 

1- 7 1-8 

2 - 6 2-6 

3- 8 3-7 

4- 5 4-5 




2. Problem 2 



5 6 7 8 



1 

2 

3 

4 



8 7 3 4 
5 5 7 4 
8 8 2 1 

9 5 9 4 



Initial Assignment Optimal Assignment 



1-7 


1-7 


2-5 


2-5 


3-8 


3-8 


4-6 


4-6 



3. Problem 3 



1 

2 

3 

4 



5 6 7 8 



5 8 4 1 

6 8 9 2 
4 2 9 1 
6 5 5 3 



Initial Assignment Optimal Assignment 



1-8 


1-8 


2-5 


2-5 


3-6 


3-6 


4-7 


4-7 



4. Problem 4 



6 7 8 9 10 



Initial Assignment Optimal Assignment 



1 

2 

3 

4 

5 



5 2 4 7 1 
3 5 9 9 4 

2 5 6 1 4 

3 9 9 3 9 
7 9 16 3 



1-10 


1-7 


2-6 


2-10 
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INTEGER X(50,50) ,AR0W(R30) , L,ACnL t VCOL,K,N, 
t MIN tPR OW , CBAR , CM! N, PC OL (230) ♦«( 230,230) ,0 
C THE NEXT 3^ LIMES R C AD TM t HF n&TA,ZSRO INVERSE &ND 

C DE t ERMINE the INITIAL ASSIGNMENT, the D A t A ARE: 

C NPROB=NUMRFR of DqnpLEMF, M=NUMR C R O <= A SS I GNE P NDnfc;, 

C X< I,J)=COS T MATRIX. NOTE t HIS PROGRAM IS DIMENSIONED TO 
C HANDLE UP TO A 15 X 15 PROBLEM. 

READ (5,10) NPROB 

C SET TIMER TO MEASURE EXECUTION TIMES. 

85 CALL SETIME 

DO 9! I =1,50 

91 A ( I ) = 0 

C A ( I ) STORES THE ASSIGNMENT FOR PRINTING PURPOSES AND TO 
C LOAD ACOL. 

DO Q2 T = 1 , ?30 
DO 92 J = 1,230 

92 R ( I , J ) = 0 
READ ( 5,10) N 
DO 1 J = 1 , N 

1 READ (5,20) ( X ( I , J ) , I =1 , N) 

WRITE (6,30) 

CALL WRTMTX ( X , N » N ) 

C WRITES INPUT nfiTA. 

TCOUNT = o 

C ICOUNT COUNTS NUMBER of ASSIGNMENTS MADE. 

2 MIN = 9999 
DO A J = 1 ,N 
DO 5 K = 1 ,N 

I c ( J.EQ.A(K)-N) GO TO A 

5 CONTINUE 

DO 3 I = 1 ,N 
I F ( A ( I ) .GT.O) GO TQ 3 
I F ( X( I , J) .GE .MIN) GO TO 3 
MIN = X( I , J) 

TROW = I 
JCOL = J 

3 CONTINUE 
A CONTINU c 

A ( TROW) = JCOL + N 
ICOUNT = ICOIJNT + 1 
I c ( ICOUNT. L T .N) GO TO 2 
WRITE (6,^0) 

DO 6 I = l ,N 

6 WRITE (6,50) I , A ( T ) 

C THF next AA LINES LOAD THE INVERT OF the PARIS (B). 

C COMMENT IS made when EACH ELEMENT OF B IS LOADED. 

I TFR = N— 1 

C I T c R COUNTS NUMBER OF ITERATIONS. 

D = N*N 
L = N*N+! 

&COL = N*N+R 
VCOL = N*N+3 
K = 1 

c load diagonal. 

DO 8 I = 1 ,n 

8 R ( I , T ) = ) 

C LOAD LAST POW EXCEPT THF CONSTANT (SIMPLEX MULTIPLIERS). 
DO 11 J = 1 , N 
DO 11 I = 1 , N 
T f ( A ( J ) - N.EO. I ) GO TO 9 
B( L » K ) = 0 

GO to 11 

9 R(L,K) = -1 
11 K = K + 1 
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C LOAD COLUMN OF CONSTANTS EXCEPT LAST ROW. 

K = l ... 

J = l 

DO 12 I = 1 ,0 
B( 1 ,L) = X ( K , J ) 

J = J + l .. . ; 

IF ( J .LE . N) GO to 12 • . • • 

K = K+1 
J = 1 

12 CONTINUE 

C LOAD ASSIGNMENT COLUMN (ACOL). 

DO 14 K = 1 ,N 
• J = (K*N) -N 
DO 13 I = 1,N 
I F ( A ( K ) - N.NE.T) GO tq 13 
B( J + T , ACOL ) = 1 

GO TO 14 .• ' 

13 CONTTNU c 

14 CONTINUE '• . 

C LOAD CONSTANT *=0R LAST ROW. ' 

DO 16 I = 1,0 

IF( B( I , ACOL) .EO.l ) GO TO 15 
B( I , VC Ol_) = 2*N +1 
GO TO 16 

15 B(T,VCOL) = (2*N + I) *(-l) 

16 CONTINUE 

SUM =0 
DO IP I = 1,0 

I R ( B( I , ACOL ) .EO. 1 ) GO T 0 17 
GO TO IB 

17 SUM = SUM + B( I ,L) 

18 CONT I N U c 

B( L , L ) = RUM * (_i ) 

C THE N c XT 23 LIMES M A K F THE R T R ST N STMdlEX ITFRATIONS. 
C THE P IR ST ° LINES FIND TH C PIVOT ROW (P<?OW) BY FINDING 
C THE minimum constant VAUI c IN EACH GROUP OF N. 

M = 1 
K = 1 

KSTAR = K + ( N-l. ) . • • 

19 MIN = 9999 

DO 21 I = K, KSTAR 

IF( B( I , L) .GE .MIN) GO TO 21 

MIN = P< I ,L) 

PROW = I • . v • -s 

21 CONttmuF • 

C URDATF VC.nt. 

B( PROWtVCOL) = M 
C PERFORM itfra t ion. 

DO 81 J = 1 , L 

DO 22 I = K, KSTAR 

I F ( I .FO. PROW) GO TO 22 

R( I , J) = B( I , J) - B( PROW, J) 

GO TO 23 

2? B< I , J) = BU , J) 

23 CONTINUE 

B( L , J ) = B ( L , J ) + B ( PROW , J ) 

81 CONTINUE 
C UPDATE indexes. 

M = M+l 
K = K + N 

KSTAR = K + (N-l) 

IF(K.LT.N*N) GO T 0 19 

C THE NEXT 24 LINES FIND T HE MINIMUM' VALU C IN the 
C TNFFASI8ILTY form FOUA t ION (CMIN). JCOMT IS A COLUMN 
C COUNTFR. JSTAR IS AN INDEX. 

49 CMIN =1 

JCONT =1 . 

JSTAR = D + 2*N 
DO 29 J = 1, JSTAR 
CBAR =0 

C NEXT 4 LINES DETFdmtnE THE APTthmftic TO BE DONE 
C DEPENDING UPON LOCATION IN STANDARD FORMAT. 
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IF(J.LE.N) GO Tn ?4 
T F ( J . L E . 2*N) GO TO ?ft 
IF(B(JCON T ,ACOL) .FO.l) GO TO 26 
CBAR = CBAR + B(L,JC0NT) 

GO TO 27 

24 K = (J-1)*N 
DO FI I = 1 , N 

51 CBAR = CBAR + B(L,I+K) 

no to 28 

25 K = J - (N+l) 

DO 53 T = 1,0, N 

5? CBAR = C RA R + B ( L , T + K ) * ( - 1 ) 

GO tp 2 ft 

26 C R AR = CBAR +R(L,JCONT) *(-l) 

27 JCOM T = JCPMT + 1 

28 IF( CBAR .GF .CMI N) GO TO 29 
CM IN = CBAR 

JCOL = J 
2 ° CONTINUE 

C THE NEXT 21 L I N r 5> TEST Tp DETFP MI NE IF AN ITERATION IS 
C RFOUTREO OR IF A CHANGE IN ASSIGNMENT IS REQUIRED. I«= NOT. 
C THE SOLUTION IS WRITTEN AND THE PROBLEM COUN T EP IS TESTFO. 
ITER = ITFR +1 
IFICMIN.LT. 0) GO TO 21 
I c ( B< L , L) .NE.O) GO tq 43 
WRITE (6,60) ITER 
K = 1 

DO B 6 I = 1,G,N 
DO R7 J = 1,N 

IF ( B(T+J-1,AC0L) .EO.0) GO to e7 
A ( K ) = N+J 
K = K + 1 
87 CONTINUE 
86 CONTINUE 

DO 8 ft I = 1 , N 
8 ft WRItc (6,50) I, AH) 

C MEASURE TIMP FROM DETECTION O c NEGATIVE CIRCUIT UNTIL 
C SOLUTION OF prpblfm. 

CALL RETIME ( I f T ) 

EL = IET * .000026 
WRITF (6,120) NRROB 
WRITE (6,130) EL 
NPROB = NRROP - l 
IF (NPROB. GT.O) GO TO 85 
STO D 

C the NEXT 2° LIMES CONSTRUCT THE COLUMN IN WHICH CMIN 
C APPEARS ( P COL ) AND c INOS T HE PIVOT ROW (PROW). 

31 I F ( JCOL.LE . N ) GO TO 34 
IFC JC0L.LF.2*N) GO TO 35 
K = JCOL - 2*N 
I C (B(K,ACPL) .EO.l) GO t 0 32 
DO 3T I = 1,0 

3T POOL ( I ) = 8(1 , K) 

GO TO 36 

32 DO 33 I = 1,0 

33 PCOL (T ) = R( I , K ) * ( -1 ) 

GO to 36 

34 K = (JCOL-1) *N 
DO 3» I = 1,0 
PCOL (I ) = 0 

DO 3ft J = 1 , N 

3ft PCOL ( I ) = PCOLU) + R (I , K+ J ) 

GO T 0 36 

35 K = JCOL -(N+l) 

DO 30 1 = i,o 
PCOL ( I) = 0 

DO 3 0 J = 1,0, N 

30 PCOLU) = PCOLU) - R(I,K + J) 

36 MIN = 9999 

00 M I = 1,0 
I F ( PCOL ( I ) .FO.l ) GO T 0 42 
GO TO 41 
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42 
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73 

74 

75 

76 



- B( PROW, J) 



IF( B( I , L) .GE .MIN) GO TO 41 
MIN = R( I ,L) 

PROW = I 

the NEXT 17 LINFS PERFORM THE SIMPLEX ITERATION. 

41 CONTINUE 

B( PROW, VCOL ) = JCOL 
DO 76 J = 1 »L 
DO 75 I = 1,0 
- IF ( I .EQ. PROW) GO to 72 

tf( pcoim .Eo.-n go to 73 

IF ( PCOL ( I ) . FO. 0 ) GO TO 72 
IF ( POOL ( I ) .FO. 1 ) GO T 0 74 . 

GO TO 75 
B ( I , J ) = B(I,J) 

GO TO 75 

B(I,J) = B ( I » J ) + B ( PROW, J ) 

GO TO 75 
B( I t J) = B( I , J) 

CONTINUE 

B(L,J) = 8 ( L * J ) - (CMIN*B( PRQW,J>) 

CONTINUE 
GO T 0 49 

r twc N c XT 25 LINES MAKFS THE CHANGE IN ASSTGMFN T . 

C GET I mf MEASURES TIME TO DETECT NEGATIVE CIRCUIT. 

43 CALL G E T I M E ( I FT ) 

EL = IET * .000026 • • 

WRITE (6,110) EL 
C RESET TIMER. 

CALL SETIMP 
DO 48 I = 1,0 

IF(B( 1 , VCOL) .GF.O) GO to 48 
DO 47 J = 1,0 

IE(B( J,ACOL).E0.1) GO TO 44 

C AROW R ?J) ^STORES ’ CONTENTS OF LAST N*N COLUMNS OF 
C ARTIFICIAL VECTOR. 

GO TO 45 

AROW(J) = B(I*J)*(-l) 

IF( AROW( J) .NE.O) GO to 46 
GO to 47 

I F ( B ( J , AC OD.EO.l) GO TO 83 
I c (B(J,ACOL) .E 0 . 0 ) GO TO 84 
GO to 47 
B( J , ACOL ) =0 
GO TO 47 
B ( J , AC OL ) = 1 
CONTINUE 

WR ITE( 6 , 80) B(I,VCOL) 

WRITE (6,90) (APOW(K) t K = 1,D) 

CONTINUE 

WR TTE ( 6 , 70 ) ITER 
GO TO 49 
FORMAT( 1 10) 

FORM AT ( 1 OT 8 ) 

30 FORMAT ( 1H1 ,10X ,' INPUT DATA') 

40 FORMAT (////) 

FORMAT (1H , I4.T6 t i- t ,T7 T4) 

FORMAT (1H0,5X’ 'PROBLEM FEASIBLE AFTER I TFRAT ION- 
114 , 5 X, ' OPTIMUM ASSIGNMENT 1^:') 

FORMAT ( 1H0,5X, ' PROGRAM INFEASIBLE AFTFR 
1 ITERATION- ' , I 4) 

FORMAT (1H0,5X, 'ARTIFICIAL VEC T OR FOP -«,T4) 
FORMAT ( 1H0,5X ,^0T 3) 

100 C QP M AT (!HO,FX, 'COS T IS:',TM _ 

110 FnR^AT (1H0,5X, 'TIME TO DETECT NEGATIVE CIRCUIT 
1 WA R ' , f I 6 . 6 , ' S c C ' ) 

120 FORMAT (1H0,5X, 'TIME TO SOLVE PROBLEM- ', 14 , 5X , 

1 ' FOLLOWS ' ) 

l^o FORMAT ( 1H0,5X,E16.6, 'SEC' ) 

END 
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SUBROUTINE WRTMTX (K,TPQW, jcou 
DIMENSION K ( 50 i 50 ) 

WRITE (6,30) (Ii I = 1»JC0L) 
PORNAT ( // 1H ,^-Xt ’ 1 / J = ' t50 T 3) 
WQITc (6,31) 

FORMAT (1H ) _ 

cQpMAT (1H , 3X , I 2 , 3X , *> 01 3 ) 

pn 33 T = 1, TROW 

WRITE (6,32) I» (K(I,J)i J = It 

RETURN 
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